Average Fidelity in n-Qubit systems 
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This letter generalizes the expression for the average fidelity of single qubits, as found by Bowdrey 
et al. [lj|, to the case of n qubits. We use a simple algebraic approach with basis elements for the 
density-matrix expansion expressed as Kronecker products of n Pauli spin matrices. An explicit 
integration over initial states is avoided by invoking the invariance of the state average under unitary 
transformations of the initial density matrix. The results have applications to measurements of 
quantum information, for example in ion-trap and NMR experiments. 



PACS numbers: 03.67.-a; 82.56.-b; 03.65.Fd 
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The fidelity, usually defined 0] as 



(the square of the Uhlmann Q formula), measures the 
correlation between two states represented by density 
matrices p and p' . The evaluation of fidelity is important 
in the storage, processing and communication of quan- 
tum information^ for example in optical qubits[||, 
ion-trap measurements or NMR experiments, espe- 
cially those with high-purity state preparation[8]. Ex- 
pression {T]) is simplified if one of the density matrices, 
say p, is in a pure state: p = |^) (\&| . Then p is idem- 
potent so that p 1 / 2 = p and 



F( P ,p') = (tiWm W7W) <*l) =<*Ip'I*> 

= ti(pp'). (2) 

This is a convenient form for calculating the fidelity of 
two output states coming from the same initial pure state 
Po = l^o) C^ol; where p is the result of a unitary trans- 
formation that preserves the purity of the state, and p' 
comes from a more general trace-preserving linear map 
M: 

p = Up rf 
P=M(p ). 



and dil — s'm8d8d(f) is an element of solid angle. Bow- 
drey et al. [l[ showed how the integration can be exactly 
h\ replaced by a finite sum in the case of a single qubit. For 
the more general n-qubit case, the corresponding integra- 
tion is difficult to envision. Nevertheless, a generalization 
to "qudit" systems of d = N independent states was de- 
veloped by Bagan et al. by expanding po hi the 
generators of SU (N) and integrating the fidelity over 
an invariant Haar measure for unitary operators. In a 
different approach exploiting global symmetries to relate 
the integral over pure states to invariant traces, Pedersen 
et al. [lfj were able to express the fidelity for any com- 
pletely positive trace-preserving linear map with a known 
Kraus-operator form , giving a formula equivalent to 
one that Emerson et al. [l2| found. In this letter, we 
present a simple approach with an explicit polarization 
basis that avoids the integration over the space of pure 
states, and with this approach we extend the Bowdrey 
et al. result to the case of n-qubit systems and show 
how a basis of pure states can be constructed for use in 
practical applications. 

To explain our approach, we first rederive the result for 
a single-qubit systems. The density matrix p of a single 
qubit can be written [l3| in terms of the polarization 
vector P, represented with components P 3 on the Pauli 
matrices 

(3) P = Yl pj °i = pVl + p2(72 + p3fT3 ' ( ? ) 



The average fidelity is defined by averaging over all pos- 
sible initial states po : 



(F) = {F(p,p')) P0 = (tr(pp')) f 



(4) 



In the case of a single qubit, the average can be found 
by integrating over the Bloch sphere [l[ 



(•5) 



where po is expanded in Pauli spin matrices o~j as 
p a = -(1 + o\ sin cos + (i2 sin 8 sin (f> + er 3 cos 8) (6) 



as 



P= 2(! + P )' 



(8) 



When p is in a pure state, P lies on the surface of the 
unit sphere 



(9) 



whereas for a general mixed state, it lies within the 
sphere: P 2 < 1. When P = 0, the state is unpolarized 
(contains no polarization information). 
If we use expression ([S]) for the initial state po, the 
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output states ([3]) are 



p=-(l + E7Ptft) 
P ' = i(M(l)+X(P)), 



(10) 
(11) 



and the fidelity becomes 

F = ~ tr [(1 + f7PC/ + )(M(l) + A4(P) 



(12) 



1 

= - tr 
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l + ^P^CTjC/t ^X(l)+^P fc X( ( 7 fc )^ 



We want to average F over initial states, that is over all 
possible directions of P. Since the coefficients P- 7 are in- 
dependent and as likely to be positive as negative, the 
average values of P- 7 and of P- 7 P fc , j ^ k, vanish. Fur- 
thermore, since po is assumed to be in a pure state, 



((P 1 ) 2 + (P 2 ) 2 + (P 3 ) 2 ) = 1, 



(13) 



and by symmetry the contribution of each component 
must be the same, 



It follows that 

{P J P k ) = ^S jk 
The average fidelity (fl"2"|) is thus 



(14) 



(15) 



The Pauli matrices can be written in terms of the pure 
states p±j = h (1 ± (Tj) as 



(Tj 1 + <7j 1 

T ~~ 2 2 
1 1 - o. 
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Pj-P 
p-p-j, 



(17) 
(18) 



where p = 1/2 is the unpolarized (maximally mixed) 
state. Then, (F) (| 16[) becomes 



< F > = \ E tr ( \ + U Pj rtM{pi) - U Pj U^M( 



P) 



(19) 



j 



The average of these expressions gives the form found by 
Bowdrey et al. [l| 

< F > = I E Tr{1 - U (P* + P-i)U f M(p)+ 

j 

Up 3 U^M{ P] ) + Up-jU^Mip-j)) (20) 

= \ E Tr<JJ Pj rtM{pj) + Up-jrfMip-j)). 
j 

This result is extended to n-qubits by defining an ap- 
propriate basis for the density matrix expansion. First 
however, we look again at the symmetry argument. 
While hardly anyone would question the equality of the 
average square components of the polarization vector in 
Eq. (fT4]) . we want to invoke a rigorous proof that does 
not depend explicitly on the integration over all states, 
and that we can extend to systems of two or more qubits. 
The proof relies on the recognition that an average over 
all initial states is invariant under unitary transformation 
of po- Such transformations are a one-to-one mappings of 
each state of the system onto another (or the same) state. 
For a single qubit, unitary transformations include rota- 
tions on the Bloch sphere [3] that could, for example, 
rotate P = PVi + P 2 a 2 + P 3 (T3 into 

P -> PPP f = P' = P 1 a 2 - P 2 cri +P 3 cr 3 , (21) 

where R is the unitary "rotor" R = exp (o20i7r/4) = 
exp (— i(737r/4) . Since the average of any expression over 
initial states is invariant under such a transformation, we 

must have (P l P 2 ) = (-P 2 P l ) and ((P 1 ) 2 ) = {{P 2 ) 

as argued above. As we see below, there is a richer se- 
lection of possible unitary transformations in systems of 
more qubits. 

For two qubits the density matrix can be expanded in 
a polarization basis {fj} as 
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E' 



(22) 



with basis elements defined in terms of a Kronecker prod- 
uct of Pauli matrices 



>ov, p,,v = 0,1,2,3, 



(23) 



where Co is the 2x2 identity matrix and a® (8) <Jq is 
excluded from the sum. The ij are traceless and satisfy 
the orthonormality condition 

trftfj-) = 4Sij . (24) 

Consequently, for any states pi , p 2 , 



tr(pip 2 ) = y^tr 
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(25) 
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If pi , P2 are the output states p,p' ([3j) , the fidelity ((2|) is 
F = ix(UpoU*M(po)) (26) 



1 

16 tr 



15 \ / 15 

1 + ]T wiUfjtf ] ( (1) + £ w k M (f fc ) 
i=i 



k=l 



For any pure state such as po, 



tr (pi) = tr (p ) 




(27) 



and thus 



E«T 



(28) 



More generally, for a state that may be mixed, 



f» 2 <3 



(29) 



with the equality holding only when the state is pure. 
The average of (|2"8"|) can be reduced by the symmetry that 
such expressions are invariant under any unitary transfor- 
mation of po, as discussed in the single-qubit case, above. 
There are unitary transformations that rotate any pair 
of polarization basis elements ij and ±ffc ([23]) into each 
other. For example, 

oi (g 1 -> (R (g 1) (o-i (g) 1) (i?t (g) l) = (j 2 <g 1 

(72 <g 1 -> (i? <g 1) (cr 2 ® 1) (i? f ® l) = -0-1 ® 1 



and 



1 ® (Ti -> (1 <g> i?) (1 ® o-i ) (1 <g i? 1 ) = 1 ® <T 2 , 



with i? = exp (i72(Ti7r/4) = exp (— ia^n/A) , are single- 
qubit rotations, whereas 

U\ (g) 1 — > -R12C1 (g) li?{ 2 = cr 2 <g CTi, 

with i?i2 = exp (<72cri ig (7i7r/4) , is a coupled-qubit rota- 
tion. Such coupled-qubit rotations are required to trans- 
form an unentangled state to an entangled one and vice 
versa. Of course single- and coupled-qubit rotations can 
be sequenced to yield still other unitary transformations. 
As a result of the invariance of state averages under all 

unitary transformations, the average of {w^ must be 
the same for every j, and we can reduce the sum for a 
pure state po by 

3HEK)Vi5(K) 2 >- ( 3 °) 
i=i 



More generally 



1 



K^ } = ^"(K) 2 } = -^ 



(31) 



Since there are also unitary transformations that sim- 
ply change the sign of the polarization components, the 
average po is unpolarized: (w J q^ = 0. Combining these 
relations, the average of the fidelity (|26|) reduces to 
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i=i 



(32) 



For the general n-qubit case, the basis element fj is a 
Kronecker product of n Pauli spin matrices cr M . The po- 
larization basis (excluding the identity element) contains 
TV 2 — 1 elements fj with N = 2", and the density matrix 
can be expanded as 



P : 



1 JV 2 -1 

3=1 



(33) 



The trace of the square of a pure-state density matrix is 
tr(po) = 1, and the fj satisfy 



tr (fjffc) = NS 



(34) 



Consequently. 
1 



1 



iV 2 



tr 



N 2 -l 



JV -1 



1 + E w & 1 + E w fa 



J'=l 



which implies 



k=l 
N 2 -l 



E {<)' 



JV -1 



E « 



= N — 1 . 



(35) 



Since, as in the 2-qubit case, there are unitary transfor- 
mations that interchange the polarization elements fj or 
change their signs, the average of (j35)) gives 

((wj)2) N_^± _J_ 
{[Wo) } ~ N*-l ~ 1+N> 

and the average including the cross terms is 

1 



-5 13 . 



"o^o/ - 1 T " ■ (36) 

The average of the fidelity F = tr(/7 / 9 t/ t A4(po)) is thus 
generalized from (|32[) to 



1 



1 



AT^-l 



N (N + 1)N 2 



[^fi^M(fj)] . (37) 
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This result confirms the qudit-fidelity relation derived 
independently by Bagan et al. 0] 

The remaining step is to express the basis elements fj 
in terms of a set of pure-state density operators. Such op- 
erators are primitive idempotents 13|, 14 1 in the algebra, 
and in the n-qubit case they can be expressed as prod- 
ucts of n commuting simple idempotents. In the 2-qubit 
case, we can use the simple idempotents 



(38) 



and note the following identities to express fj = <Jp,® o v 
as a linear combination of pure-state density matrices 

® <J V = (P +fll/ - P-^u) (P +UfM + P- Vtl ) , n > v (39) 
= (P +u0 - P-po) (P+Ou - P-Ou) , V > 0, v > 

A similar procedure is possible for the general n-qubit 
case. Although the explicit expressions become rather 
lengthy, the procedure is straightforward. 
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